arXiv:l502.00279vl [hep-ph] 16 Jan 2015 


The lowest hidden charmed tetraquark state from QCD sum rules 

Zhi-Gang WangQ 

Department of Physics, North China Electric Power University, Baoding 071003, P. R. China 

Abstract 

In this article, we study the SS type scalar tetraquark state cqcq in details with the QCD 
sum rules by calculating the contributions of the vacuum condensates up to dimension-10 in 
the operator product expansion, and obtain the value Mz c = (3.82^'g|) GeV, which is the 
lowest mass for the hidden charmed tetraquark states from the QCD sum rules. Furthermore, 
we calculate the hadronic coupling constants Gz c tj c it and Gz c dd with the three-point QCD 
sum rules, then study the strong decays Z c —> r/ c n , DD, and observe that the total width 
T z a ~ 21MeV. The present predictions can be confronted with the experimental data in the 
futures at the BESIII, LHCb and Belle-II. 
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1 Introduction 

The scattering amplitude for one-gluon exchange in an SU(N C ) gauge theory is proportional to 

tkitfj =- ^ (SjkSu - SikSji) + (SjkSu + 5 ik 5ji ), ( 1 ) 

where the t a is the generator of the gauge group, and the i,j and k, I are the color indexes of the two 
quarks in the incoming and outgoing channels respectively. For N c = 3, the negative sign in front 
of the antisymmetric antitriplet indicates the interaction is attractive and favors the formation of 
the diquark states in the color antitriplet, while the positive sign in front of the symmetric sextet 
indicates the interaction is repulsive and disfavors the formation of the diquark states in the color 
sextet pL. 

The antitriplet diquark states have five Dirac tensor structures, scalar C 75, pseudoscalar (7, 
vector <77^75, axial vector C 7 M and tensor Ccr MJ/ . The structures C 7^ and Ca^ u are symmetric, the 
structures C75, C and <77^,75 are antisymmetric. The attractive interactions of one-gluon exchange 
favor formation of the diquarks in color antitriplet 3 0 , flavor antitriplet 3/ and spin singlet l s (or 
flavor sextet 6/ and spin triplet 3 S ) Ell, so the favored configurations are the scalar and axial- 
vector diquark states. The scalar (S') and axial-vector (A) heavy-light diquark states have almost 
degenerate masses from the QCD sum rules [4][5j. In Refs.|6][7j, we take the C75 —C75, — 

Cjfj ,75 — C 7^75 type interpolating currents to study the masses of the scalar tetraquark states in a 
systematic way using the QCD sum rules, and observe that the SS and AA type scalar tetraquark 
states have almost degenerate masses, about 4.36 GeV, which is much larger than that from the 
phenomenological models mmm- 

In Ref.[8], Ebert, Faustov and Galkin calculate the masses of the excited heavy tetraquarks 
with hidden charm within the relativistic diquark-antidiquark picture based on the quasipotential 
approach, and obtain the values Mj—q = 3.852 GeV and 3.812 GeV for the AA and SS type scalar 
tetraquark states cqcq, respectively. While L. Maiani et al obtain the values Mj- 0 = 3.832 GeV 
and 3.723 GeV for the AA and SS type scalar tetraquark states cqcq respectively in the type-I 
diquark-antidiquark model [Rj, and Mj = 0 = 4.000GeV and 3.770GeV for the AA and SS type 
scalar tetraquark states cqcq respectively in the type-II diquark-antidiquark model |10j . In those 
model-dependent studies, the masses of the AA-type scalar tetraquark states are larger than that 
of the SS type scalar tetraquark states. 

1 E-mail: zgwang@aliyun.com. 


1 





In Refs.illj T2, TTi. [T3J Q3} SI], we explore the energy scale dependence of the hidden charmed 
(bottom) tetraquark states and molecular states in details for the first time, and suggest a formula 

M = y/M 2 x/y/z - ( 2 Mq )2 , ( 2 ) 

with the effective heavy Q-quark mass Mq to determine the energy scales of the QCD spectral 
densities in the QCD sum rules, which works well. According to the formula, the energy scale 
[i= 1 GeV taken in Refs. mm is too low to result in robust predictions. 

In Ref. [13], we choose the — C^ v type interpolating currents to study the AA-type scalar, 
axial-vector and tensor tetraquark states in details with the QCD sum rules. The predicted masses 
of the axial-vector and tensor tetraquark states favor assigning the Z c ( 4020) and -Zc(4025) as the 
J PC = 1+- or 2 ++ diquark-antidiquark type tetraquark states. While there are no experimental 
candidates to match the predicted mass of the scalar tetraquark state M J=0 = (3.85+° Jg) GeV. 
The value is consistent with the prediction Mj =0 = 3.852 GeV based on the quasipotential approach 
H], while the upper bound reaches the prediction Mj = o = 4.000 GeV based on the type-II diquark- 
antidiquark model [TO] , According to Refs. nsnni, the SS -type scalar tetraquark states have 
smaller masses than that of the corresponding AA-type scalar tetraquark states. It is interesting 
to see whether or not such conclusion survives when confronted with the QCD sum rules. In 
Refs.fTT j H2I fl3l . we observe that the masses of the SA or AS type axial-vector tetraquark states 
are larger than that of the AA type scalar tetraquark states. So the SS scalar tetraquark state 
maybe the lowest tetraquark state. 

In this article, we study the scalar SS-type hidden charmed tetraquark state (thereafter we will 
denote it as Z c ) by calculating the contributions of the vacuum condensates up to dimension-10, 
and try to obtain the lowest mass based on the QCD sum rules. Furthermore, we calculate the 
hadronic coupling constants Gz^tv and Gz c dd with the three-point QCD sum rules, then study 
the strong decays Z c —»• t] c tt , DD. 

The article is arranged as follows: we derive the QCD sum rules for the mass and pole residue 
of the scalar tetraquark state Z c and for the hadronic coupling constants Gz cT)c7r and Gz c dd in 
section 2; in section 3, we present the numerical results and discussions; section 4 is reserved for 
our conclusion. 


2 QCD sum rules for the scalar tetraquark state 

In the following, we write down the two-point correlation function II(p) in the QCD sum rules, 

n(p) = i J d 4 xe ipx {0\T{J(x)J^(0)}\0), (3) 


J(x) = e ijk e imn u j (x)C'y 5 c k (x)d m (x)'y 5 Cc n (x ), 


( 4 ) 


where the i, j , k , m, n are color indexes, the C is the charge conjugation matrix. 

At the hadronic side, we can insert a complete set of intermediate hadronic states with the 
same quantum numbers as the current operator J(x) into the correlation function II(p) to obtain 
the hadronic representation [T7i nn nn . After isolating the ground state contribution of the scalar 
tetraquark state, we get the following result, 


n(p) 


A| 

M z c - p 2 


( 5 ) 


where the pole residue Xz c is defined by (0| J(0)\Z(p)) = Xz c ■ 
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In the following, we briefly outline the operator product expansion for the correlation function 
II(p) in perturbative QCD. We contract the it, d and c quark fields in the correlation function II(p) 
with Wick theorem, and obtain the result: 


i 4 xe ipx 


n(p) = i e vk e imn e i'j'k' e i'm'n' J d 

Tr \ lb C kk ' {x)^CU^' T {x)c\ Tr \ l5 C n ' n {-x)^CD m ' mT (-x)C 


( 6 ) 


where the Uij(x), Dij(x) and Cij{x ) are the full u, d and c quark propagators respectively (the 
Uij(x) and Dij(x) can be written as Sij (a:) for simplicity), 


Sij ( x) 


i5ij Sij(qq) 

2tt 2 x 4 12 
8ijX 4 {qq){g 2 s GG) 
27648 


Sijx 2 ( qg s aGq ) W s G a a0 t^ tya 01 ? + a a/3 $ 

192 32ir 2 x 2 

- ^(djcr pv q z )a^ - ^(qj^qi) 7/* + • ■ • , 


idjjX 2 #g 2 (gq) 2 

7776 


( 7 ) 


Cij(x) = 


yAa/3 


d 4 ke~ ik ' x 


9sG™pt™j a a P(ft + rn c ) + {fi + m, c )a ai3 


(27t) 4 J “ \f6 — m c 4 (k 2 — m 2 ) 2 

g s D a G^(f^ a + f x ^) g 2 s {t a t%G a a/3 G b ^(f a0,xv + f a ^ u + f a ^) 

3(fc 2 — m 2 ) 4 4(fc 2 — to 2 ) 5 

Qt + m c ) 7 A (/(f + m c ) 7 “$ + m c )^(]t + m c ), 

(# + m c )j a Q( + m c ) 7 /3 (ft + m c ) 7 M (# + m c )Y(ft + m c ), 


( 8 ) 


and t n = 4j-, the A 71 is the Gell-Mann matrix, DC = d a — ig s G™t n ITS], then compute the integrals 
both in the coordinate and momentum spaces to obtain the correlation function II(p) therefore the 
QCD spectral density. In Eq.(7), we retain the terms {qjU^qi) and (qj^^qi) originate from the 
Fierz re-arrangement of the {qi.qj) to absorb the gluons emitted from the heavy quark lines so as 
to extract the mixed condensate and four-quark condensate ( qg s uGq ) and g 2 (qq} 2 , respectively. 

Once the analytical expression is obtained, we can take the quark-hadron duality below the 
continuum threshold so and perform Borel transform with respect to the variable P 2 = —p 2 to 
obtain the following QCD sum rule: 


where 



/4m? 


ds p(s ) exp 



( 9 ) 


p(s) = Po(s) +P 3 (s) + P4,{s) + p 5 (s) + Pe{s) + P7(s) + Ps(s) + Pio(s), (10) 


i py f /*1— y 

Po (s) = 5127r g J dy J dzyz{l-y - z) 3 (s-m 2 c ) 2 (7 s 2 - esml+mi) , (11) 

P '3 00 = ~ J d V J z dz(y + z)(l-y-z)(s-m%)(2s-ml) , ( 12 ) 


ml ,a s GG, 


rvf 


r^-y 


M = “384 dy 


dz I -3 + A ) (1 - y- s ) 3 

yZ Z Z 


2s — m 2 + -jp-8 (s — ml) 

1 rv (7(7 r y f G~ v 

[(— 7 :—)j dy J dz (y + z) (1 - y - z) 2 (lOs 2 - 12sm 2 c + 3rn^) , (13) 


5127T 4 ' 7T 
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m c {qg s <jGq) f Vf f 1 v ( 

Pb{s) = -77 —i - / dy dz (y + z) (3s - 2m c ) 


647T 4 


'Vi 


r ' dy rdz ( i + 1) (1 - „ - *) (3, - ■ml) , 

v / 


647T 4 


( 14 ) 


Pe(s) = 




12tt 2 


*+mr«r 


>Vi 

9 2 s (qq) 2 [ vs 

512tt 4 


108tt 4 




dz yz ^ 2s — m 2 + -jr-8 (s — ml) | 


ri-y 


dy 


dz (1 - y - z) \ 2 ( - + -J (3s - 2ml) + + z i 


m 2 c [2 + m 2 c S (s-m 2 )]} 

“H? f dv L, v * (1 -* - z) { 3 (; + f) (3s - 2 ^> + (? + 

m 2 [2 + m 2 c 5 (s — to 2 )] + (y + z) [12 (2s — ml) + 2 mlS (s — We)] } > 


2 2 
Z Z 


( 15 ) 


P7(s) = 


Wc(ot) ,a s GG 


288tt 2 


7T 


rvf 

) dy I 

j Vi J Zi 


x ~ v ( V z 1 1 

( — -f 3 H 2 ^ 2 ) (^ — 2 / — Z ) 

yz ' 5 ^ z ^ 5 1 


1 +pr) J ( s - m =) 

m c {qq ) ^ a s GG ^ f Vf ^ 


rl-y 


967r 2 




^ ^ ) C 1 — 2 / — ^) { 2 + (« - } 


m c {qq) ,a s GG f Vf 

-96p-<—>A ^ 


rl— 1 / 


dz {2 + mlS (s — ml )} 


m c (qq) a s GG f Vf fn 2 w ~2\i 

"5767T 2- 7T / dy l 2 + m ^( S_m Gl 


(16) 


Ps(s) = - 


+ 


m 2 c(qq)(qg s crGq) 
24t r 2 

m l(qq) {qgs&Gq) 


dy [1 + ^)5 {s- ml) 


Pw(s) = 


48tt 2 

ml(qg s oGq ) 2 f 1 
1927t 2 T 6 ./o 

ml(qq) 2 ,a s GG 


' dy— -7 

0 2/(1 27) 


d (s-m 2 ) , 


-<- 


+ 


216T 4 ' 7T 
ml(qq) 2 a s GG 


72T 2 

ml(qg s aGq) 
192t r 2 T 4 


7T 

2 rl 


dyml 5 (s — ml) 

>jf 


(1 _„)3| { (»-"*“) 


'0 ir (i-y) 2 


, m 2 c (qg s aGq) 2 f 1 ^ 

4 77T - ;—^7^— / dy 


dy , 2 

y 

1 

0 ~ w y(i-y) 

1 


d (s - TO 2 ) 


dy 


+ 


384tt 2 T 2 y 0 y(l — y) 
ml(qq) 2 a s GG 

1 I 

0 


?7l 2 d (s — TO 2 ) 
d (s - TO 2 ) 


<- 

216T 6 x 7T 


-) / dy ml d (s — to 2 ) 

./O 


(17) 


(18) 
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the subscripts 0, 3, 4, 5, 6 , 7, 8 , 10 denote the dimensions of the vacuum condensates, y/ = 

l+\/l"K/s _ l-\A-4 ml/s __ _ yml 2 _ (.V+z)m 2 c ~ 2 _ m* rVf , , fl , 

- - - m c - - - -i m c ~ „n j J Vi a V ^ Jo a y ’ 


2 > Hi — 2 > — ys-m'i ’ "‘c — ^ > ""c — y(l-y) ’ J Vi ' J 0 

f_ y dz->-f 0 y dz when the (5 functions 5 (s — m(!) and 5 (s — to))) appear. We take into account 
the vacuum condensates which are vacuum expections of the operators of the orders O(a^) with 
k < 1 consistently. 

Differentiate Eq.(9) with respect to y?, then eliminate the pole residues A z c , we obtain the 
QCD sum rule for the mass of the scalar tetraquark state, 


Zi = 


~ O 

mi = 


M l = 


f4m 2 c ds d(-l/T 2 )P( S ) ex P ( 72 ) 

IZ 2 ds p ( s ) ex p (-#) 


(19) 


In the following, we perform Fierz re-arrangement to the current J both in the color and 
Dirac-spinor spaces to obtain the result, 


J = 


—ccdu + ciy$c diy^u — cy^cdy^u — cy 11 y$c dy ^y^u + —ca^cda^u 


+cu dc — ciyyu diy^c + C'y^u dy^c + cy^y^u dy^y^c — -ca^u da^c 


( 20 ) 


the components couple to the meson pairs XcO^q (980), r] c n + , J/ipp + , Xci7r + , XciZ (1260), h c hf (1170), 
(D o (2400)S o (2400))+, (DD)+, {D*D*)+, (D 1 (2420)£>i(2420)) + , (D 1 (2430)Si(2430))+, respec¬ 
tively. The strong decays 

Zf{ 0 ++ ) -t Xco4 (980), 7 ?c tt ± , J/^p ± , X ci^ ± , Xciof (1260), h c hf (1170), (D o (2400)D o (2400)) ± , 
(DD)* , (D*5*) ± , (D 1 (2420)Di(2420)) ± , (D 1 (2430)D 1 (2430)) ± , ( 21 ) 


are Okubo-Zweig-Iizuka super-allowed, if they are kinematically allowed. The diquark-antidiquark 
type tetraquark state can be taken as a special superposition of a series of meson-meson pairs, and 
embodies the net effects. The decays to its components (meson-meson pairs) are Okubo-Zweig- 
Iizuka super-allowed, but the re-arrangements in the color-space are non-trivial [201121| . 

The numerical analysis indicates that the ground state mass of the SS -type scalar tetraquark 
state is about 3.82 GeV, the strong decays 

Zt( 0 ++ ) -l JfcTT* , Xd^ , {DD)* , (22) 

are kinematically allowed. The decay Z±( 0 ++ ) —> Xd'> T± takes place through relative P-wave and 
is kinematically suppressed. 

Now we write down the three-point correlation functions III (p, q ) and II 2 (p, q ) to study the 
strong decays Z±(0 ++ ) —¥ rj c 7 r ± , ( DD ) ± , 

ni(p, q) = i 2 J d i xd i ye ip ' x e iq " y (0\T { J I?o (a;) J w (p) J(0)} |0), 
n 2 (p,g) = i 2 Jd 4 xd 4 ye ip x e iq y (0\T{J D -(x)J D o(y)J(0)}\0), (23) 

where the currents 

J Vc (x) = c.(x)iyoc{x ), 

Jn{y) = u(y)iy 5 d(y), 

J D -{x) = c(x)iy 5 d(x ), 

JD°(y) = u(y)iy 5 c(y), (24) 
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Figure 1: The connected Feynman diagram contributes to the correlation function Hi(p, q), where 
the dashed and solid lines denote the heavy quark and light quark lines, respectively. Other 
diagrams obtained by interchanging of the heavy quark lines or light quark lines are implied. 



Figure 2: The connected Feynman diagram contributes to the correlation function II 2 (p, q ), where 
the dashed and solid lines denote the heavy quark and light quark lines, respectively. Other 
diagrams obtained by interchanging of the heavy quark lines and (or) light quark lines are implied. 


interpolate the mesons ?y c , 7 r, D ~, D°, respectively. 

We insert a complete set of intermediate hadronic states with the same quantum numbers as 
the current operators into the three-point correlation functions H\(p,q) and H- 2 (p,q) and isolate 
the ground state contributions to obtain the following results, 

fnM%fr) c M^ c \z c Gz c ri c n —q • p 

2 (m u + m d )m c ( M z c ~ P' 2 )(Kc - P 2 )( M % - q 2 ) 

f1,MfjXz c Gz c DD _ -g-p _/ 2 rx 

(m c + m q ) 2 {M 2 o -p/ 2 ){Ml-^){M 2 D -q 2 ) + "' ’ 1 J 


IIi (p,q) = 
n 2{p,q) = 


where p' = p+q , the fo, fr] c and /, r are the decay constants of the mesons D , p c and w, respectively, 
the Gz cVc ir and Gz c dd are the hadronic coupling constants. In the following, we write down the 
definitions, 


<0|J, C (0)1^0?)) 
(0|J w (0)|7r(g)> 

(0\J D (0)\D(p/q)) 

{vc(p)^(q)\Z c {p')) 

(D(p)D(q)\Z c (p')) 


UK 

2 m c ’ 

UK 

m u +m d 
fpMl ^ 
m c + m q 

iq • pGz c r) c n{q ) 5 
-iq ■ pG Zc DD(q 2 ) ■ 


(26) 

(27) 


We carry out the operator product expansion and take into account the color connected Feyn¬ 
man diagrams mm, and obtain the following results, 


Hi (p,q) 


m c (qg s aGq) f 1 q-p 

32ir 2 q 2 J 0 X m 2 — x[l — x)p 2 


(28) 
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112 (p, g) 


1 


m c (qg s aGq) q ■ p 
64w 2 q 2 — m 2 
m c (qg s aGq) q ■ p 
647r 2 m 2 — p 2 


dx 

dx 


1 + x 

m 2 — (1 — x)p 2 
2 — x 
xq 2 — m 2 


(29) 


In Fig.l and Fig.2, we draw the connected Feynman diagrams contribute to the correlation func¬ 
tions IIi(p, q) and II 2 (p, g), respectively. The H±(p,q) and n 2 (p, q) can be expanded in terms of 

the cos 9, n l / 2 (p ; q) = II 0 (p 2 , q 2 ) + H 1 (p 2 , q 2 ) cos 9 + II 2 (p 2 , q 2 ) cos 2 9 -1 -, at the QCD side, where 

the 9 is the included angle of the Euclidean momenta p and q, i.e. cos 9 = p ■ q/y/q 2 p 2 . There 
exists only one term (II 1 (j> 2 , q 2 ) cos9) for the IIi(p,q), while there exist two terms (II 0 (p 2 ,g 2 ) and 
H 1 (p 2 ,q 2 )cos9) for the n 2 (p, q). At the phenomenological side, the hadronic coupling constants 
Gspp'(p, q) have the possible forms Gg PP ,, Gg PP , cos 9, G 2 SPP , cos 2 9 , • • •, where the S denotes 
the scalar mesons, the P and P' denote the pseudoscalar mesons. In the present case, it is bet¬ 
ter to choose the form Gg PP , cos 9, as the correlation functions II 1 (p, g) and II 2 (p, q) both have 
the term proportional to cos0 at the QCD side. The cos 9 is the pertinent tensor structure, as 
the correlation functions III (p,g) and n 2 (p, q) should have the same tensor structure at the phe¬ 
nomenological side. There exists some shortcoming, if we choose the form G° SPP , and take the 
replacement 2p ■ q = p' 2 — p 2 — q 2 , then set p 2 = p' 2 and perform the Borel transform with respect 
to the variable P 2 = —p 2 , as the p , q and p' are not independent variables, the cos 9 cannot be 
replaced. 

Once the analytical expressions of the correlation functions Hi (p,q) and II 2 {p,q) at both the 
QCD side and hadron side are obtained, we perform the Borel transform with respect to the variable 
P 2 = —p 2 by setting p 2 = p' 2 , then take the quark-hadron duality and obtain the following QCD 
sum rules, 


UM 2 f Vc M 2 c X Zc G Zc , c7I \ ( M 2 c \ f M 2 Zc \\ 

2 (m u + m d )m c (M 2 Zc M 2 J { 6XP { T 2 ) &XP { T 2 )) 

m c {qg s aGq) Q 2 + M 2 f 1 1 / m 2 \ 

327 r 2 Q 2 J 0 x(l — x) \ x(l — x)T 2 ) 


G ZcVc7 t exp 



(30) 


fpM%\z c G Zc DD 
{m c + m q ) 2 {M 2 -M 2 ) 


exp 


Ml 

2^2 


exp 




2 1 2 


+ CZcDD exp (—^) 


{q2 + Ml) 


647T 2 
2 — x 


dx 


I + 2 : 


Q 2 + m 2 (1 - x) 


exp 


(1 -x)T 2 


xQ- 


■ m% 


■ exp 


2 1 2 


(31) 


where the so is the continuum threshold parameter for the Z c , and the C Zcric ^ and C Zc dd are 
unknown parameters introduced to take into account the single-pole contributions associated with 
pole-continuum transitions. In numerical analysis, we will denote the right sides of Eqs. (30-31) as 
F\(Q 2 ) and F 2 (Q 2 ) respectively. In the three-point QCD sum rules, the single-pole contributions 
are not suppressed if a single Borel transform is taken. 


3 Numerical results and discussions 


The vacuum condensates are taken to be the standard values (qq) = —(0.24±0.01 GeV) 3 , ( qg s aGq) = 
m o(qq), m o = (0-8 ± 0.1) GeV 2 , ( a ^ G ) = (0.33 GeV) 4 at the energy scale n = 1 GeV 17, .18,19, 
1221123] . The quark condensate and mixed quark condensate evolve with the renormalization group 


equation, ( qq){p ) 


(<w)(Q) 


«„(Q) 1 9 
«s(m) 


and ( qg s aGq)(p) 


{qg s aGq)(Q) 


a s (Q) 

a s (p) 


2 

27 
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The hadronic input parameters are taken as M n = 0.13957 GeV, /„■ = 0.130 GeV, M D ± = 
1.8695 GeV, M D o = 1.86491 GeV, f D = 0.208 GeV, M Vc = 2.9837 GeV, f Vc = 0.350 GeV [24j[25l 

[26] . 

We take the values m u (fi = lGeV) = rrid(n = lGeV) = m q (n = lGeV) = 0.006GeV from the 
Gell-Mann-Oakes-Renner relation, and choose the MS mass m c (?n c ) = (1.275 ± 0.025) GeV from 
the Particle Data Group [23], and take into account the energy-scale dependence of the MS masses 
from the renormalization group equation, 


m q (n) 

m c (n) 

Mm) 


m q ( lGeV) 


s(v) 


2 {m c ) 


a s (lGeV) 
12 


_a s (m c )_ 


1 

b Q t 


1 _ id log t bl(\og 2 t- log t- 1) + b 0 b 2 


b Q t 




(32) 


where t = log£, b 0 = 332 " f , 


^ _ 153 — 19 n/ 

0l — 24 if 5 : 


62 — 


2857 - 5 ®^ + ^! 

128 tt 3 


A = 213 MeV, 296 MeV 


A 2 ’ 12 tt 

and 339 MeV for the flavors n/ =5, 4 and 3, respectively 

Now we study the mass and pole residue of the SS type scalar tetraquark state. We impose the 
two criteria (pole dominance and convergence of the operator product expansion) on the hidden 
charmed tetraquark state to choose the Borel parameter T 2 and threshold parameter so- 

In the heavy quark limit, the c (and b) quark can be taken as a static well potential, which binds 
the light quark q' to form a diquark in the color antitriplet channel or binds the light antiquark 
q to form a meson in the color singlet channel (or a meson-like state in the color octet channel). 
Then the heavy tetraquark states are characterized by the effective heavy quark masses Mq (or 


constituent quark masses) and the virtuality V = M^/y/z; — (2Mq) 2 (or bound energy not as 

robust). It is natural to take the energy scale n = V, the formula works well for the V(3872), 
Z c (3885), Z c { 3900), Z c (4020), Z c { 4025), Z(4050), Z{ 4250), V(4360), £(4430), V(4630), F(4660), 
Zb( 10610) and Z&( 10650) in the scenario of tetraquark states [TT1 IT21H31 HI] . The relation 


/Y/Z = (2M C ) 2 + fj, 2 , (33) 

with the value M c = 1.8 GeV determined in previous works mmmm puts a strong constraint 
on the masses of the possible tetraquark states. 

The mass gaps between the ground states and the first radial excited states are usually taken 
as (0.4 — 0.6) GeV, for example, the £(4430) is tentatively assigned as the first radial excitation of 
the £ c (3900) according to the analogous decays, 

£ c (3900) ± —► J/tjj^ , 

£(4430 ) ± -> (34) 

and the mass differences M Z ( 443 o) — ^z c ( 3900 ) = 576 MeV, M., 1,1 — Afj/p = 589 MeV flOl [27[ [28]. 
The relation 


y/so — Mx/y/z + (0-4 — 0.6) GeV, 


(35) 


puts another strong constraint on the masses of the possible tetraquark states. 
In calculations, we observe that 


/i t M z i , 
M 4- M z t, 


( 36 ) 

















T 2 (GeV 2 ) 


Figure 3: The pole contribution with variations of the Borel parameter T 2 and threshold parameter 
so, where the A, B , C, D , E, F denote the threshold parameters y'so = 4.0, 4.1, 4.2, 4.3, 4.4, 
4.5 GeV, respectively. 


from the QCD sum rule in Eq.(19). While Eq.(33) indicates that 

A* t M z t, 

At 4- Mz l ■ (37) 

There must be a compromise, which leads to the optimal energy scale p, mass Mz and threshold 
parameter sq. 

In Fig.3, the contribution of the pole term is plotted with variations of the threshold parameter 
so and Borel parameter T 2 at the energy scale /i = 1.3 GeV. From the figure, we can see that the 
value y^So <4.1 GeV is too small to satisfy the pole dominance condition and result in reasonable 
Borel window. 

In Fig.4, the contributions of different terms in the operator product expansion are plotted with 
variations of the Borel parameter T 2 for the threshold parameter yGo = 4.3 GeV at the energy 
scale p = 1.3 GeV. From the figure, we can see that the Do , Do, Do, Do and Dg, where the D; 
denote the contributions of the vacuum condensates of dimensions D = i, play an important role, 
while the D 4, D 7 and Dio play a minor important role. At the value T 2 < 2.0 GeV, the D3, D 5 , 
Do and Dg decrease monotonously and quickly with increase of the T 2 , which cannot lead to stable 
QCD sum rules. At the value T 2 = (2.2 — 2.6) GeV 2 , D 3 |D 5 | ^ Do^> |D 8 | and Di 0 -C 1%, 
the operator product expansion is well convergent, although D 0 ~ 20%. We approximate the 
continuum spectral density by Pqcd(s)Q(s — Sq)', the contributions of the quark condensate (qq) 
and mixed condensate ( qg s aGq ) can be very large. 

In this article, we take the Borel parameter T 2 = (2.2 — 2.6) GeV 2 , the continuum threshold 
parameter ^/£o = (4.2 — 4.4) GeV and the energy scale p = 1.3 GeV, the pole dominance is well 
satisfied. The Borel parameter, continuum threshold parameter and the pole contribution are 
shown explicitly in Table 1. The two criteria (pole dominance and convergence of the operator 
product expansion) of the QCD sum rules are fully satisfied, furthermore, the relations in Eq.(33) 
and Eq.(35) are also satisfied. 

Taking into account all uncertainties of the input parameters, finally we obtain the values of 
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Figure 4: The contributions of different terms in the operator product expansion with variations 
of the Borel parameter T 2 , where the 0, 3, 4, 5, 6, 7, 8, 10 denote the dimensions of the vacuum 
condensates. 


J pa 

T 2 (GeV 2 ) 

^(GeV) 

pole 

M z ( GeV) 

Az 

0 ++ 

2 .2-2.6 

4.3 ±0.1 

(49 - 74)% 

q QQ+U.U8 

O.OZ—o 08 

!.79+^4 x 10- 2 GeV D 


Table 1: The Borel parameter, continuum threshold parameter, pole contribution, mass and pole 
residue of the scalar tetraquark state. 



T 2 (GeV 2 ) 


Figure 5: The mass with variations of the Borel parameter T 2 . 
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Figure 6: The pole residue with variations of the Borel parameter T 2 . 


the mass and pole residue of the SS type scalar tetraquark state, which are shown explicitly in 
Figs.5-6 and Table 1. 

The central value of the present prediction Mz c = (3.82^q'q|) GeV for the SS type scalar 
tetraquark state is smaller than that of the AA type scalar tetraquark state Mj = o = (3.85+ 0 09 ) GeV 
obtained in Ref. m- The predictions based on the QCD sum rules are consistent with the val¬ 
ues Mj— o = 3.852 GeV and 3.812 GeV for the AA and SS type scalar tetraquark states cqcq 
respectively from the quasipotential approach [ 8 ]. 

Now we take the mass Mz c and pole residue A z c as basic input parameters to study the 
hadronic coupling constants Gz c ri C Tt and Gz c dd , and take the same threshold parameter and Borel 
parameter as in the QCD sum rule for the mass and pole residue. In calculations, we choose the 
unknown parameters as Gz c r / r tt = 0.0009 GeV and Cz c dd — 0.0004 GeV 6 to obtain stable QCD 
sum rules with variations of the Borel parameter T 2 at the Borel windows T 2 = (2.2 — 2.6) GeV 2 ; 
the left side and right side of the QCD sum rules coincide. In fact, it is not necessary to choose the 
same Borel parameters both in the two-point and three-point QCD sum rules. If we take larger 
Borel parameter, say T 2 = (2.5 — 3.0) GeV 2 instead of T 2 = (2.2 — 2.6) GeV 2 , we should alter the 
unknown parameters C'^ eI)c7r and Cz c dd slightly, then obtain stable QCD sum rules, the resulting 
values of the hadronic coupling constants change slightly. 

Based on Eqs.(30-31), we can study the Q 2 dependence of the right side of the QCD sum rules, 

Fi(Q 2 ) oc y wl| (38) 

at the region of large (or intermediate) Q 2 due to the tiny mass of the ^ r, while the ^(Q 2 ) has no 
such simple Q 2 dependence due to the heavy quark mass m c and heavy meson mass Mjj. In the 
limit Q 2 —» oo, 


f 2 (q 2 ) 


m c (qg s aGq) f 1 

——Jo dX 


(1+x 
\ 1 — X 


exp 


(1 - z)T 5 


2 — x 

H-exp 

x 



, (39) 


which is independence on Q 2 . In Fig.7, we plot the central values of the ^(Q 2 ) with variations 
of the Q 2 at the range Q 2 = (1 — 5) GeV 2 for the Borel parameters T 2 = 2.2 GeV 2 , 2.4 GeV 2 
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Figure 7: The central values of the .^(Q 2 ) with variations of the Q 2 . 



Figure 8: The central values of the hadronic coupling constants with variations of the Q 2 , where 
the A and B denote the Gz c dd(Q 2 ) and Gz cr)c ^{Q 2 ), respectively. 
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and 2.6 GeV 2 , respectively. From the figure, we can see that the Q 2 dependence of the F 2 (Q 2 ) is 
rather mild and can be neglected approximately. The left sides of the QCD sum rules in Eqs.(30- 
31) have no explicit Q 2 dependence, the Q 2 dependence is embodied in the right sides of the 
QCD sum rules ( F\{Q 2 ) and F 2 (Q 2 )), so the hadronic coupling constants Gz c r) c ir and Gz c dd are 
independent on the Q 2 in the limit Q 2 —> oo, the conclusion survives even for much smaller Q 2 , 
say Q 2 = (1 — 5) GeV 2 according to Eq.(38) and Fig.7. The central values of the Gz cVc tt(Q 2 ) and 
Gz c dd{Q 2 ) can be fitted to the following constant forms, 

G Zc r,AQ 2 ) = 0.43 GeV" 1 , 

G Zc dd(Q 2 ) = 1.06 GeV -1 , (40) 

at the region Q 2 = (1 — 5) GeV 2 ; the uncertainties of the Gz c rj c n and Gz c dd are about 25% and 
18%, respectively. We plot the central values of the hadronic coupling constants Gz c dd(Q 2 ) and 
Gz c n c n{Q 2 ) with variations of the Q 2 at the region Q 2 = (1 — 5) GeV 2 for the Borel parameter 

T 2 = 2.4 GeV 2 in Fig.8. From the figure, we can see that the fitted functions in Eq.(40) are 

satisfactory. We extend the coupling constants to the physical regions without difficulty, and 
calculate the partial decay widths, 


?c 7T 

Y z c -+DD 


where 


Vr} c TT 


PDD 


c 2 

^ Z C T) C -K 


(M 2 Zc - Ml 


- Mlf 


Pr/cK 


32ttM f c 

(M 2 c - M 2 d+ - M 2 d0 ) 2 Vdd 


C 12 

^Z C DD 


32?rMf 


(3.0 ± 1.5) MeV , 

= (17.9 ±6.4) MeV, 


\J \Ml - (M Vc + M„) 2 \ [Ml - (M Vc - M n y 
2 AD 


yJ\Ml - (M d+ + M d o) 2 ] [Ml - (M d+ - M d o) 

2M~Zc 


(41) 


(42) 


The total width Tz c of the Z c ( 3820) can be approximated by Tz c ^ Vc -k + T z c ^dd, the numerical 
value is about (20.9 ± 6.6) MeV. The radiative decay widths can be estimated by assuming vector 
meson dominance, for example, T z ±^ p ± oc a\Y z ±^j^, p ± | for the radiative decays Z^( 3820) —> 

J/ip*p± —f "yp ± , the partial decay widths are of the order O(KeV) due to the factor a = . 

The strong decays 2^(3820) —> J/ipp± are kinematically forbidden, the values of the T z ±^j/^* p ± 
are complex, so we take \Y z ±^j,^, p± \. The contributions of the radiative decays to the total width 
Fz c are small and can be neglected. 


4 Conclusion 

In this article, we calculate the contributions of the vacuum condensates up to dimension-10 in 
the operator product expansion, study the SS type scalar tetraquark state cqcq in details with 
the QCD sum rules. In calculations, we search for the optimal Borel parameter and threshold 
parameter to satisfy the energy scale formula M z = (2M C ) 2 ± p 2 and the experiential threshold 
formula y/so = Mz + (0.4 — 0.6) GeV, where the p is the energy scale of the QCD spectral density, 
and obtain the values Mz c = (3.82toos) GeV and Xz c = (1.79+Q24) x 10 _2 GeV 5 . The central 
value of the mass of the SS type scalar tetraquark state is smaller than that of the AA type scalar 
tetraquark state, the SS type scalar tetraquark state cqcq maybe the lowest hidden charmed 
tetraquark state. Furthermore, we calculate the hadronic coupling constants Gz c ti c -k and Gz c dd 
with the three-point QCD sum rules by taking into account the color-connected diagrams, then 
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study the strong decays Z c —> i] c ir , DD , and observe that the total width Fz c ~ 21MeV. The 
present predictions can be confronted with the experimental data in the futures at the BESIII, 
LHCb and Belle-II. 
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